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Symmetries in neutrino physics are explored using analogies to fermion pairing in many-body 
systems. In particular, the SO(5) symmetry of the most general neutrino mass Hamiltonian with 
both Dirac and Majorana mass terms as well as the invariants of certain limits of the collective 
neutrino Hamiltonian are discussed. The latter Hamiltonian finds applications in the Early Universe 
and core-collapse supernovae. 
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Recent years have seen an increasingly intense activity in neutrino physics. Not only has neutrino mixing been 
experimentally established, but also the three mixing angles have been measured with reasonable accuracy. Direct 
£NJ \ neutrino mass measurements have not yet reached down to the neutrino mass scale chosen by the Nature, but 
differences between squares of masses are deduced from oscillation experiments again with reasonable accuracy despite 
the fact that the mass hierarchy remains an open question. 

Symmetries play a crucial role is constraining the behavior of physical systems. In this contribution I describe 
how symmetry concepts apply to neutrino physics using analogies to fermion pairing in many-body systems. One use 
of this analogy, namely the SO(5) symmetry which connects the Dirac and Majorana masses of neutrinos with the 
see-saw mechanism is discussed in the next section. 

Since our current knowledge of neutrino parameters is rather robust, we can now address applications of neutrino 
physics in astrophysical settings with more confidence. One such application is the self-interacting neutrino gas present 
in the core-collapse supernovae and the Early Universe. Symmetries, the integrability and associated constants of 
motion of the Hamiltonian describing such a many-neutrino gas is described in the section after the next one. 
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II. SO(5) ALGEBRA AND THE NEUTRINO MASS 



Both the Dirac mass term 

= m D J <Px(4> L ip R + h.c.) (1) 
as well as the left- and right-handed Majorana mass terms 

H% =H^ + H* = \m L J d 3 x(^ L + h.c) + l -m R J d 3 x(M R + h.c), (2) 

where the charge-conjugate spinor is defined as tp = Cip T , can be written in terms of the generators of an SO (5) 
algebra [1]. The generators of the SU{2)d subalgebra of this SO(5) are 



£>_ = / d'xiML) = D\, D = - J d 6 x(tPlij L - ^r), (3) 



whereas the generators of the SU(2)l X SU(2)r subalgebra are 



and 



L+ = \f d 3 x(i> L i>i) = Ll,L Q = - A j d 3 x(^ L - M[) (4) 



R+ = \f d 3 x(r R yb R ) =Rl,R = l I d 3 x(M R - ^y f: i ( r, ) 
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Note that Dq is not an independent operator, but the sum of Lq and Rq. With the addition of two more operators, 



A+ = / d A x \-r L C l0 ^ R \ ,A_= d A x 



f R loC(^ L ) T \ , (6) 

one has the ten generators of the 50(5) pairing algebra. 

To explore the physical meaning of this S0(5) algebra, one can consider the Pauli-Gursey transformation 

i> -> V' = + bj 5 ^ c , |a| 2 + |6| 2 = l, (7) 

which describes particle-antiparticle mixing. The two operators A + and A- of Eq. (j6|) along with Aq = Rq — Lq 
generate another SU(2) algebra that we call SU(2)pq. The most general element of the associated SU(2)pq group is 

fj — e -r* A- e -log(l + \ T \ 2 )A & tA + e iipA /g-j 

under which the field ip transforms as 

i > -^// = U^= J^L^ -t*^ c \. (9) 
Clearly this is a Pauli-Gursey transformation with 

e *v/2 _ T * p iip/2 

a ■■ 



The most general neutrino mass Hamiltonian: 

H m = m D (D+ + DJ) + m L (L + + L_) + m R (R+ + 

sits in the SO(5)/SU(2)pg X U(l) x coset where U(l) x is generated by Do = Lq + Rq. Under the Pauli-Gursey 
transformation Dirac mass terms would transform into a mixture of Dirac and Majorana mass terms, providing an 
algebraic framework for the see-saw mechanism. 

Although the discussion here pertains only to one single flavor, it is also possible to generalize these algebraic 
arguments to multiple neutrino flavors @. 

III. SYMMETRIES AND INTEGRABILITY OF THE SELF-INTERACTING NEUTRINO GAS 

We first consider only two flavors of neutrinos: electron neutrino, v e , and another flavor, v x . Introducing the 
creation and annihilation operators for a neutrino with three momentum p, we can write down the generators of the 
neutrino flavor isospin algebras Q: 

J+(p) = 4(pH(p)> j -(p) = <4(pK(p)> 

Mp) = i(4(pK(p)-a f e (pH(p)). (10) 

The integrals of these operators over all possible values of momenta generate the global flavor isospin algebra. Using 
the operators in Eq. (|10[) the Hamiltonian for a neutrino propagating through matter takes the form 



H u = J d3 P^- cos20J o (p) + isin20(J + (p) + J_(p)) - V2G F J d 3 p N e J (p) 



(11) 



In Eq. dill) , the first integral represents the neutrino mixing and the second integral represents the neutrino forward 
scattering off the background matter. In writing this equation a term proportional to identity is omitted as such 
terms do not contribute to the neutrino oscillations. Neutrino-neutrino interactions (see references Q through (2TJ] 
are described by the Hamiltonian 

2T„„ = V2^ y"d 3 pd 3 q(l-costf pq ) J(p) ■ J(q), (12) 

where $ pq is the angle between neutrino momenta p and q and V is the normalization volume. Note that the presence 
of the (1 — cosi? P q) term in the integral above is crucial to recover the effects of the Standard Model weak interaction 
physics in the most general situation 1 . Note that in the extremely idealized case of isotropic neutrino distribution 



1 For a recent discussion of the impact of the physics beyond the Standard Model on matter-enhanced neutrino oscillations see [2lj | and 
references therein. 
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and a very large number of neutrinos, this term may average to a constant and the neutrino-neutrino interaction 
Hamiltonian simply reduces to the Casimir operator of the global SU(2) algebra. This is unlikely for core-collapse 
supernovae, but is the case for the Early Universe [22l l23l|. 

The discussion above pertains to a gas comprised of neutrinos only. Inclusion of antineutrinos as well requires 
introduction of a second set of SU(2) algebras. Similarly for three flavors two sets of SU(3) algebras are needed [24j . 
Both extensions are straightforward, but tedious. 

Defining the auxiliary vector quantity 

B = (sin 26, 0, - cos 20), (13) 
the total Hamiltonian with two flavors, containing one- and two-body interaction terms, can be written as 

fftota. = H v + H vv = ( J2 ^T-£ • J p - V2G F N e J Q p ) + ^ (1 - cos # pq ) J p • J q (14) 

V p / p,q 

The evolution operator for the system represented by the Hamiltonian in Eq. (|14[) 



i^r = (Hu + H vv ) U, (15) 
at 

can be approximately evaluated using the stationary phase approximation to its path integral representation 4] . This 
is equivalent to reducing the Hamiltonian H ulJ , to a one-body one in an RPA-like approximation. In this approximation 
the product of two commuting operators 0\ and O2 is approximated as 

6 1 6 2 ~ Oi(o 2 > + (Oi)o 2 - (6 1 ){6 2 ), (16) 

where the expectation values should be calculated with respect to a well-chosen state \ty) which satisfies the condition 
(0102) = (Oi)(C>2) ■ One then obtains the single-angle Hamiltonian 

H ~ H RPA = u pB -Jp + P-J- (17) 

p 

In writing Eq. (fl~7|) matter terms are neglected and the polarization vector P was defined as 

P P , S = 2<J p>s ). (18) 

Using the SU(2) coherent states associated with the flavor isospin in calculating the operator averages in the above 
equations yields the standard reduced collective neutrino Hamiltonian, widely used in the literature [4|. 

Both the full Hamiltonian and its one-body reduction possess an SU{N) f rotation symmetry in the neutrino flavor 
space [U, HH, [26[. Such a complex nonlinear system could presumably exhibit further symmetries. A few of the 
conserved quantities in collective neutrino oscillations were already noted in the literature [3, H] ■ To look for further 
invariants, defining \i — , r = fit, and uj p — ^^j- [13], one can write the Hamiltonian of Eq. (TI"4"|) in the 

single-angle approximation as 

H = ^uj p B ■ J p + J '■ J. (19) 

p 

In writing Eq. (|19[) , it was assumed that neutrino- neutrino interaction term is dominant and hence the matter terms 
are ignored. This Hamiltonian preserves the length of each spin 



Lp — Jp • Jp , 



H, Lp 



, (20) 
as well as the total spin component in the direction of the "external magnetic field", B 



Co = B ■ J 



H,C 



. (21) 



The conservation law depicted in Eq. (|20[) is nothing but the conservation of the total number of neutrinos with a 
given momentum as neither neutrino mixing nor coherent forward scattering of neutrinos off one another change that 
number. 
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The Hamiltonian in Eq. (fT9|) is similar to the reduced BCS Hamiltonian of the many-body theory. (However, one 
should note that the sign of the pairing term is opposite to that in the BCS Hamiltonian). One can exploit this 
duality to uncover its symmetries. Eigenstates of the reduced BCS Hamiltonian were written by Richardson using a 



Bethe ansatz in [28| and later generalized by Gaudin [29, 30|. Since the BCS Hamiltonian considered by Richardson is 
integrable, there are constants of motion associated with it (3l| . Using this analogy one can write down the constants 
of motion of the collective neutrino Hamiltonian in Eq. (fT4| as [32] 



L=B- L + 2 Y Jp ' Jq . (22) 
z — ' uj„ — UJ n 

p q 

The individual neutrino spin- lengths discussed above, L p , are independent invariants. They are set by the initial 
conditions and are not changed by the evolution of the system under the collective Hamiltonian. However one has 
Cq = Y] p hp. The Hamiltonian itself is also a linear combination of these invariants: 



H = w p hp + y]L p 



Including antineutrinos, the conserved quantities for each neutrino energy mode p take the form 

^=B-J P + 2V^A + 2V^-, (23) 

where for antineutrinos we defined cup = — Conserved quantities hp for different antineutrino energy modes are 

ft p 4j p + 2yM + 2 yM. (24 ) 

The invariants of the one-body Hamiltonian of Eq. (|17|) can be written from those as 

I p = 2(k p ) =B-P p + V iWjL + ^iWA (25 ) 

and 

I p = 2(hp)=B-Pp+ V EPlI± + yEl^L. (26 ) 

It was shown that existence of such invariants could lead to collective neutrino oscillations [33[ • 

As in the BCS theory, the one-body Hamiltonian of Eq. (| 1 T[) does not conserve particle number. Particle number 
conservation can be enforced by introducing a Lagrange multiplier: 

ff RPA -> H nrA + ljJ° . (27) 

Diagonalization of this Hamiltonian gives rise to the phenomena called spectral split or swapping in the neutrino 
energy spectra [13, [3^ - l36j with the Lagrange multiplier playing the role of the swap frequency. 

The fact that invariants of the full Hamiltonian are also invariants of the one-body Hamiltonian of Eq. p7|) when 
they are properly linearized provides confidence in the aptness of the linearization procedure itself. One should also 
note that another linearization procedure has been used to carry out flavor-stability analysis of dense neutrino streams 



IV. CONCLUSIONS 



Neutrinos play an important role in astrophysical settings, including core-collapse supernovae [38|. Core-collapse 
supernovae are one of the possible sites for the r-process nucleosynthesis. Yields of r-process nucleosynthesis are 



5 



determined by the electron fraction, or equivalently by the neutron-to-proton ratio, n/p. Interactions of the neutrinos 
and antineutrinos streaming out of the core both with nucleons and seed nuclei determine the n/p ratio. Hence it is 
crucial to understand neutrino properties, neutrino interactions, and symmetries of those interactions in a many-body 
environment. 

We examined the many-neutrino gas both from the exact many-body perspective and from the point of view of 
an effective one-body description formulated with the application of the RPA method. In the limit of the single 
angle approximation, we showed that both the many-body and the RPA pictures possess many constants of motion 
manifesting the existence of associated dynamical symmetries in the system. The existence of such constants of motion 
offer practical ways of extracting information even from exceedingly complex systems. Even when the symmetries 
which guarantee their existence is broken, they usually provide a convenient set of variables which behave in a relatively 
simple manner depending on how drastic the symmetry breaking factor is. The existence of such invariants naturally 
lead to associated collective modes in neutrino oscillations. It should be emphasized that existence of invariants does 
not obviate numerical analysis as the stability of such collective behavior still needs to be numerically studied. 

Whether there are invariants associated with the more realistic multi-angle collective neutrino Hamiltonian is an 
open question. Note that even the multi-angle picture may need to be modified because of the presence of neutrinos 
that undergo direction-changing scattering outside of the neutrinosphere . 
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